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Abstract. In this note, we give a generalization of the inversion formulas of 
Pestov-Uhlmann for the geodesic ray transform of functions and vector fields on 
simple 2-dimensional manifolds of constant curvature. The inversion formulas 
given here hold for 2-dimensional simple manifolds whose curvatures close to 
a constant. 



1. Introduction 

Let {M, dM, g) be a C°° compact Riemannian manifolds with boundary. A 
variant of the classical Radon transform on Euclidean space is the geodesic ray 
transform on Riemannian manifolds defined as follows: 

/m/(7)= / {.f{n{t)),T{t))dt, 



where 7 : [0, Z(7)] — > Af is a maximal geodesic parameterized by arc length and m 
indicates the rank of the symmetric tensor field / G L'^{M). 

We will be interested only in the cases m = (functions) and m ~ 1 (vector 
fields) in this note and denote their geodesic ray transforms by /q and Ii respec- 
tively. 

The geodesic ray transform is not injective in general. One needs additional 
restrictions on the metric and one such restriction is to assume that the Riemannian 
manifold (M, dM, g) is simple [Sha94 defined as follows: 

Definition 1. A compact Riemannian manifold with boundary is simple if 

(a) The boundary dM is strictly convex: (V^i^, ^} < for ^ G Tx{dM) where v 
is the unit inward normal to the boundary. 

(b) The map exp^ : exp~^ M — * M is a diffeomorphism for each x G M . 

It is known that on a simple Riemannian manifold, /q/ uniquely determines / 
and Iif uniquely determines the solenoidal component of /. For references to these 
works, we refer the book of Sharafutdinov |Sha94| . Then, similar to the classical 
Radon inversion formula, it is natural to ask whether there exists explicit inversion 
formulas for a function or a vector field in terms of its geodesic ray transform. In 
general this is a hard problem and such formulas are known in only in special cases 
plilOQ). 

In |PU04j , Pestov and Uhlmann found Fredholm-type inversion formulas for the 
geodesic ray transform of functions and vector fields for simple 2-dimensional mani- 
folds. These formulas become exact inversion formulas for 2-dimensional manifolds 
of constant curvature, even when conjugate points are present along geodesies. 
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A brief remark regarding notation. We use notation that is standard in integral 
geometry literature. Ours is consistent for the most part with |Sha94| . In this note, 
SM is the unit sphere bundle and r(x, ^) is the length of the maximal geodesic 
starting at x £ dM in the direction ^ e d+SM := {(x,^) e dSM : {vix),£,} > 0}. 

The Fredholm-type inversion formulas of Pestov-Uhlmann are given by the fol- 
lowing theorem: 

Theorem 1. |PU04|, Theorem 5.4] Let {M,g) be a 2- dimensional Riemannian 
manifold. Then 

f + ^'f-^ ^^±1*1 {o'*H{Iof)-\o^sM) , / e L\M). 

h + {W*fh = {a*H{hH±h) + \o^sM) , h e H^(M). 

Here W is the operator (W* is its adjoint) on L'^{AI) defined by 



with 

As shown in |PU04j . for manifolds of constant curvature, W = W* = and 
hence these formulas becomes exact inversion formulas. 

In this note, we generalize these formulas to simple 2-dimensional manifolds 
whose curvatures are close to a constant. We show that in this case, the inversion 
formulas are given by convergent Neumann series expansions. Our main result is a 
generalization of the above result: 

Theorem 2. There exists aC > such that if M is a simple 2-dimensional man- 
ifold with Gaussian curvature K such that \\VK\\cq < C, the following inversion 
formulas hold: 

f = (/ + (^5±/r (a*i?(/o/)"|a+SM)) , ./ e L\M). 

h={I+ (^l-I* {a*H{hn±h)+\o^sM)^ , h e Hl(M). 

Remark: The proof relies on getting bounds for the operator W (and hence W*) 
in terms of the gradient of the curvature K. Hence we will not give definitions 
of the terms appearing on the right hand side of the formulas in Theorems [1] or [5] 
which can be found in Pestov-Uhlmann's papers |PU05[ IPU04j . 

As shown in |PU04j . W is a smoothing integral operator extendible as a map 
yy : L^{M) -> C°°(M) with kernel, 



Wi^x,y) ^ -Q(x,exp,^(,)) ''^--^(-P-^^;(^^^'v^ . (1) 

The function Q (equations ^ and (O) and the partial differential operator de 
(equation Q) are defined in the appendix. 
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2. The proof 

We prove the foUowing lemma. Here the derivatives are with respect to time. 
The functions a and b are defined in the appendix; see equation 

Lemma 1. Let (p — bdga — adgb. Then denoting — K o ^ , if satisfies the 
following ordinary differential equation, 

(^(3) + AK^if' + 2K'^if = ~2dgK^, 

with initial conditions, 

(p(0)-^'(0)-^"(0)-0. 

Proof. First of all we have 

ab' -a'b=l. (2) 
For, let (p = ab' — ba'. Then (p' = ab" — a"b. From equation ([5]) we get that (p' — 
and so (/) is a constant. Since ((>{{)) = 1, we have the claim. With this we now show 
that If — bdga — adgb satisfies the ODE above. 

if' = b'dea + bdea' — a'dgb — adgb' . 

From ^ we get 

b'dga + adb' - bdga' ~ a'dgb = 0. 

This gives 

if' = 2{b'dga - a'dgb) = 2{bdga' - adgb'). 
Differentiating again, we get 

If" = 2{b"dga + b'dga' - a"dgb ~ a'dgb'). 

Using equation this reduces to 

if" = 2{-K^if + b'dga' ~ a'dgb'), 

where K-y = K o ^. Differentiating yet again, and as in the steps above, we finally 
get, 

if^^^ + 2Ky'f + AK^if' = -2dgK^, (3) 
ft now follows directly from these equations that if{Q) — "^'(0) = </?"(0) = 0. □ 

Notation: The norm jj • || in the proof below denotes the sup norm unless indicated 
otherwise. Also in order to avoid proliferation of subscripts, we will use the same 
letter C to denote different constants. 

Proof of Theorem [H We now prove the main theorem. We rewrite equation ([3]) as 
a first order differential equation. We get 

1 0\ / 

1 U2 + 





( 









\-2k: 




where 

ipi ^ if, if 2^ ip' and if-i = f". 
For simplicity, let us write this as a system of the form 

X'{t)^A{t)X{t)+B{t), 
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where X,B and matrix A depend also on {x,S, j). From [Cod61j . since ^(0) — 0, 
we have a solution of this differential equation to be 

t 

X{t) = / ^~\s)B{s)ds. 



where $ is the fundamental matrix of the homogeneous differential equation, 

X'{t) ^ A{t)X{t). 

Since the manifold is compact, we have ||$||, < oo. From the relation, 

deK^ = {C±,VK^), 
and using the fact that SM is compact, we have a C such that 

\deK^\<C\\VK\\. 
Combining these inequalities we get, 

\^'{x,tt)\ < \X{x,^,t)\<Ct\\VK\\. 

for some C > 0. Since 

t 

fit) = J 'P'is)ds, 



we have 

Mt)\<Ct^VK\\. 

We can initially work with (p"{t) and by the same arguments as above, we also get 

W)\<Ct^S7Kl 

for a different constant C. 

Since the manifold is simple, we have b for t 0, since b{0) = 0. Now we 
write b{t, x, y) = tb{x, y, t) with 6 7^ 0. Therefore for a suitable C > 0, 

\q{x,^,t)\ = \de^\<Ct\\VKl 

where the norm of y K is the sup norm. Since tQ{x, t^) — q{x, ^, i), we have 

\Q{x,tO\<C\\VK\\. 

Since the remaining terms in 



u./ N r)t -1/ AAdet(exp^i)'(a;,y)y^(^ 
W{x,y) = -Q(a;,exp^ (y)) ^— - 

V5(y) 

are bounded above by compactness of Af , we have 

||P^||<C||VX||. 

Therefore we have 

l|W|U2_c~(M) <C||Vif||. 

So now choosing HV-fCH to be small enough, we have |1>V|| < 1. Hence we have 
inversion formulas involving Neumann series expansions recovering the function 
from its geodesic ray transform. A similar argument works for the recovery of the 
solenoidal part of a vector field from its geodesic ray transform. This completes the 
proof of the theorem. □ 
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Appendix A. The kernel of W 

For completeness and because the function q defined in equation ([6]) is critical 
for the proof of Theorem [21 we sketch below, Pestov-Uhlmann's |PU04] derivation 
of the integral kernel of the operator W. 

Recall that the operator W is defined as 

where 



2n 



For a function u on SM, is defined as 



We have 



= ^("°^^)l|«l=i' where = (a;,C/|C|). 

SxA/ 

We define two Jacobi vector fields along the geodesic 'y{x,£,,t) as follows: Let 
x(s), — e < s < e be a curve starting at x in the direction Now parallel translate 
the vector ^ along this curve, call it ^(s) and consider the variation by geodesies, 
j{x{s),^{s),t). The vector field 

is a Jacobi vector field along 7 with the following initial conditions, 

X{x,tO)^^^, DtX{x,tO) = 0. 
It can be also be written as 

x{x,^,t) = n±jix,^,t). 

We now define another Jacobi vector field by considering the variation by geodesies, 
7(0;, ^(s), i), where f(s) is a smooth curve in SxM with initial tangent vector ^j^. 
The Jacobi vector field 

dg-f{x,tt) ■■= Y{x,tt) = -^\s=Qlix,(is),t) (4) 

as 

has initial conditions, 

y(x,e,0)=0, DtY{x,C,0)^U- 

Since X and Y are vector fields normal to 7 and because of dimensional reasons 
these two fields must be proportional to the parallel translate of the vector ^± along 
the geodesic 7. Let this parallel translate be denoted 7^. Then there exists two 
smooth functions a(x,^,t) and b{x,^,t) such that 

X — a7j_, Y — 67 j_. 

The functions a and b satisfy the scalar Jacobi equations, 

a" + Ka = b" + Kb = 0, (5) 
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with 

a{x,tO) = l, a'(x,^,0)=0, 5(x,^,0) = 0, b'ix,^,0)^l. 
We now write, 



r(x..a 
S^M 

S^M 



We now define a function q on 

G = {{x, ^, t) : (a:, C) G -t(x, -C) < * < ^G^, 0, ^ 0} 

by 

q{x,U)^de[^). (6) 

Using this we define a function Q E C°°{TAI) by, 

Q{x,t0^tq{x,^,t). (7) 

The existence of this function Q follows from the fact that the geodesic jix, ^, t) is 
smooth as a function of {x, tS^) and the initial conditions of the Lemma [T] We now 
get the integral kernel W{x, y) in equation ([T]) by a change of variables involving 
the inverse of the exponential map. 
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